An Infrared Re normalization Group Limit Cycle in QCD 
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Small increases in the up and down quark masses of QCD would tune the theory to the critical 
renormalization group trajectory for an infrared limit cycle in the three-nucleon system. At critical 
values of the quark masses, the deuteron binding energy goes to zero and the triton has infinitely 
many excited states with an accumulation point at the 3-nucleon threshold. The ratio of the binding 
energies of successive states approaches a universal constant that is close to 515. The proximity of 
physical QCD to the critical trajectory for this limit cycle explains the success of an effective field 
theory of nucleons with contact interactions only in describing the low-energy 3-nucleon system. 
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The development of the renormalization group (RG) 
has had a profound effect on many branches of physics. 
Its successes range from explaining the universality of 
critical phenomena in condensed matter physics to the 
nonperturbative formulation of quantum field theories 
that describe elementary particles [lj . Most simple appli- 
cations of the renormalization group involve fixed points 
that are approached either in the infrared limit or in the 
ultraviolet limit. However, the renormalization group can 
be reduced to a set of differential equations that define 
a flow in the space of coupling constants. A fixed point 
where the flow vanishes is only the simplest topologi- 
cal feature that can be exhibited by such a flow. An- 
other topological feature is a limit cycle, which is a 1- 
dimensional orbit that is closed under the RG flow. The 
possibility of RG flow to a limit cycle was proposed by 
Wilson in 1971 2J. Glazek and Wilson have recently con- 
structed a simple discrete Hamiltonian system that ex- 
hibits RG flow to a limit cycle 3J. However, few physical 
applications of RG limit cycles have emerged. 

The purpose of this letter is to point out that Quantum 
Chromodynamics (QCD) is close to the critical trajectory 
for an infrared RG limit cycle in the 3-nucleon sector. It 
can be tuned to the critical trajectory by small changes 
in the up and down quark masses. The proximity of the 
physical quark masses to these critical values explains the 
success of a program initiated by Efimov to describe the 
3-nucleon problem in terms of zero-range forces between 
nucleons |£| . An effective-field-theory formulation of this 
program by Bedaque, Hammer, and van Kolck exhibits 
an ultraviolet RG limit cycle |5| . The proximity of phys- 
ical QCD to the critical trajectory for an infrared limit 
cycle shows that the ultraviolet limit cycle of Ref. || is 
not just an artifact of their model. 

In the late 1960's, Wilson used the renormalization 
group (RG) to explain universality in critical phenom- 



short-distance degrees of freedom while leaving the long- 
distances physics invariant define a RG flow on the multi- 
dimensional space of coupling constants g for all possible 
terms in the Hamiltonian: 
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ena [1|. Transformations of a system that integrate out 



where A is an ultraviolet momentum cutoff. Standard 
critical phenomena are associated with infrared fixed 
points g* of the RG flow, which satisfy /3(g*) = 0. The 
tuning of macroscopic variables to reach a critical point 
corresponds to the tuning of the coupling constants g to 
a critical trajectory that flows to the fixed point g* in 
the infrared limit A — ► 0. One of the signatures of an 
RG fixed point is scale invariance: symmetry with re- 
spect to the coordinate transformation r — * Ar for any 
positive number A. This symmetry implies that dimen- 
sionless variables scale as powers of the momentum scale, 
perhaps with anomalous dimensions. 

In a classic 1971 paper, Wilson suggested that the 
renormalization group might also be relevant to the 
strong interactions of elementary particle physics Q ■ At 
that time, the fundamental theory for the strong interac- 
tions had not yet been discovered. It was believed to 
involve quarks, and hints that the strong interactions 
might have scaling behavior at high energies had been ob- 
served in experiments on deeply-inelastic lepton-nucleon 
scattering. Wilson suggested that simple high-energy be- 
havior can be explained by simple RG flow of the rele- 
vant coupling constants in the ultraviolet limit A — > oo. 
The simplest possibility is RG flow to an ultraviolet fixed 
point. Another simple possibility is RG flow to an ultra- 
violet limit cycle. A limit cycle is a 1-parameter family 
of coupling constants g*(0) that is closed under the RG 
flow and can be parametrized by an angle < 9 < 2ir. 
The RG flow carries the system around a complete orbit 
of the limit cycle every time the ultraviolet cutoff A in- 
creases by some factor Ao- One of the signatures of an RG 
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limit cycle is discrete scale invariance: symmetry with re- 
spect to the coordinate transformation r — > AJJr only for 
integer values of n. This symmetry implies that dimen- 
sionless variables with vanishing anomalous dimensions 
are periodic functions of the logarithm of the momentum 
scale with period ln(Ao). The fundamental field theory 
for the strong interactions, QCD, was eventually discov- 
ered. QCD has a single coupling constant a 8 (A) with 
an asymptotically-free ultraviolet fixed point: a s (A) — > 
as A — > 00 6J. The renormalization of QCD does not 
involve an ultraviolet limit cycle. 

We now turn to an independent theoretical develop- 
ment whose connection with RG limit cycles was first 
pointed out in Ref. Q. In 1970, Efimov discovered a re- 
markable effect that can occur in the 3-body sector for 
nonrelativistic particles with a resonant short-range 2- 
body interaction || . If the particles are identical bosons 
and if the resonance is in the S-wave channel, the strength 
of the interaction is governed by the S-wave scatter- 
ing length a. Efimov showed that if \a\ is much larger 
than the range rg of the interaction, there are shallow 
3-body bound states whose number increases logarithmi- 
cally with |a|/ro- In the resonant limit a — ► ±00, there 
are infinitely many shallow 3-body bound states with an 
accumulation point at the 3-body scattering threshold. 
The ratio of the binding energies of successive states 
rapidly approaches the universal constant A 2 , ~ 515, 
where Ao = e*l s ° « 22.7 and sq ~ 1.00624 is a transcen- 
dental number. In subsequent work, Efimov studied both 
the bound state spectrum and low-energy 3-body scatter- 
ing observables for identical bosons with large scattering 
length [9j. He showed that the observables for different 
values of a are related by a discrete scaling transforma- 
tion in which a — > Ag a, where n is an integer, and lengths 
and energies are scaled by the appropriate powers of Ag . 

The Efimov effect can also occur for fermions with 
equal masses and large scattering lengths if they have 
at least 3 distinct spin states. Nucleons are examples of 
fermions with large scattering lengths. The spin-singlet 
and spin-triplet np scattering lengths are a s = —23.8 fm 
and at = 5.4 fm. They are both significantly larger than 
the NN effective range, which is r = 1.8 fm in the spin- 
triplet channel. Efimov used this observation as the basis 
for a qualitative approach to the 3-nucleon problem 0. 
His starting point was to take nucleons as point particles 
with zero-range potentials whose strengths are adjusted 
to reproduce the scattering lengths a s and at- The ef- 
fective range and higher order terms in the low-energy 
expansions of the phase shifts were treated as perturba- 
tions. This approach works well in the 2-nucleon system, 
giving accurate predictions for the deuteron binding en- 
ergy. This is no surprise; it simply reflects the well-known 
success of the effective range expansion in the 2-nucleon 
system Remarkably, Efimov's program also works 

well in the 3-nucleon system at momenta small compared 
to m^. In the triton channel, the Efimov effect makes it 



necessary to impose a boundary condition on the wave- 
function at short distances. The boundary condition can 
be fixed by using either the neutron-deuteron scattering 
length or the triton binding energy as input. But if one of 
them is used as input, the other is predicted with surpris- 
ing accuracy. The accuracy can be improved by taking 
into account the effective range as a first-order perturba- 
tion [n|. 

The Efimov effect was revisited by Bedaque, Hammer, 
and van Kolck within the framework of effective field the- 
ory (EFT) |l2( . The problem of bosons with mass m and 
large scattering length a can be described by a nonrel- 
ativistic field theory with a complex-valued field ip and 
Hamiltonian density 

H = {n 2 /2m)Vip* ■ Vip + <? 2 (A)0/;*V) 2 . (2) 

For convenience, we set h = 1 in the following. In the 2- 
body sector, the nonperturbative solution of the field the- 
ory can be obtained analytically. Renormalization can be 
implemented by adjusting the 2-body coupling constant 
32(A) as a function of the ultraviolet cutoff A so that 
the scattering length is a. Other 2-body observables are 
then independent of A and have the appropriate values 
for bosons with zero effective range. 

In the 3-body sector, the nonperturbative solution of 
the field theory can be obtained by solving integral equa- 
tions numerically. These integral equations have unique 
solutions only in the presence of an ultraviolet cutoff 
A. The resulting predictions for 3-body observables, al- 
though finite, depend on the ultraviolet cutoff and are 
periodic functions of ln(A) with period tt/so. Bedaque, 
Hammer, and van Kolck showed that the quantum field 
theory could be fully renormalized to remove the residual 
dependence on A in the 3-body sector by adding a 3-body 
interaction term g3(A)(tp*ip) 3 to the Hamiltonian density 
in J2J) [l2| • The dependence of 3-body observables on the 
cutoff decreases like 1/ A 2 if the 3-body coupling constant 
has the form 33(A) = -4mg 2 (A) 2 F(A)/A 2 , where the di- 
mensionless function H(A) has the form 

= sin[s ln(A/A«) - arctan(l/s )] 
^ ' sin [s In (A/ A*) + arctan(l/s )] 

for some value of A* . With this renormalization, 3-body 
observables have well-defined limits as A — > 00, but they 
depend on the parameter A* introduced by dimensional 
transmutation. Note that H(A) is a periodic function 
of ln(A). Thus the renormalization of the field theory 
involves an ultraviolet limit cycle. 

An effective field theory (EFT) of nucleons with con- 
tact interactions only has also been applied to the 3- 
nucleon problem 0, 0|. This contact EFT involves 
an isospin doublet N of Pauli fields with two indepen- 
dent 2-body contact interactions: aiN c N c ^ <7iN and 
N^T k N c N c ^T k N, where N c = o- 2 t 2 N* . Renormalization 
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FIG. 1: The inverse scattering lengths l/at and l/o s as func- 
tions of nirr as predicted by the EFT with pions of Ref . [l^l . 

in the 2-body sector requires that the two coupling con- 
stants be adjusted as a function of the ultraviolet cut- 
off A to obtain the correct spin-singlet and spin-triplet 
scattering lengths a s and at- Renormalization in the 
3-body sector requires the 3-body contact interaction 
<7{N C N cjf tjj N i\H <JiUj N with a coefficient proportional 
to Q. Thus the renormalization involves an ultraviolet 
limit cycle. 

The low-energy few-nucleon problem can also be de- 
scribed by an EFT that includes explicit pion fields as 
well as contact interactions between the nucleons. The 
renormalization of such an EFT does not involve any RG 
limit cycle. This suggests that the RG limit cycle in the 
EFT of Ref. @ is just an artifact of the model. We will 
argue to the contrary that it is actually a hint of an in- 
frared limit cycle in QCD. 

Our argument is based on recent work in which an 
EFT with explicit pions was used to extrapolate nuclear 
forces to the chiral limit of QCD flilllsT ]. In this limit, the 
masses m u and rrid of the up and down quarks are zero 
and the pion is an exactly massless Goldstone boson asso- 
ciated with spontaneous breaking of the chiral symmetry 
of QCD. According to the chiral extrapolation of Epel- 
baum, Meifiner, and Glockle ^J, the small binding en- 
ergy 2.2 MeV of the deuteron is a fortuitous consequence 
of the physical values of m u and ma- When extrapo- 
lated to the chiral limit, the deuteron has a much larger 
binding energy comparable to the scale \/{m,Nr1) ~ 10 
MeV set by the NN effective range. Conversely, if ex- 
trapolated farther from the chiral limit, the deuteron's 
binding energy decreases to and then it becomes un- 
bound. This effect is illustrated in Fig. ^ which shows 
the chiral extrapolation of the inverse scattering lengths 
l/a t and 1/eu as functions of m w from Ref. |l5j. In the 
EFT of Ref. |l5j, the coefficients of some of the 2-nucleon 
contact interactions are not well-constrained by the low- 
energy data. The bands in Fig. ^ are obtained by vary- 
ing those coefficients over natural ranges. The width of 
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FIG. 2: The binding momenta k = (mB 3 ) 1/2 of pnn bound 
states as a function of m n calculated using the contact EFT of 
Ref. 0- The circles indicate the triton ground state and ex- 
cited state. The crosses give the binding energy of the physical 
deuteron and triton, while the dashed lines give the thresh- 
olds for decay into a nucleon plus a deuteron (left curve) or a 
spin-singlet deuteron (right curve). 



the error band of course shrinks to zero at the physical 
value of tyi-k . The prediction of Ref. is that the crit- 
ical value m^t at which l/a t = is in the range 170 
MeV < m^j < 210 MeV, which is not much larger than 
the physical value of m^. The extrapolation of l/a s has 
larger uncertainties. It may increase to zero at some crit- 
ical value m„. )S greater than 150 MeV, in which case the 
spin-singlet deuteron is bound for > nru^, or l/a s 
may remain negative. Beane and Savage |l4| have ar- 
gued that the extrapolation errors in the chiral limit are 
larger than estimated in Ref. 0] . For the small extrap- 
olation to the region of larger where the deuteron 
becomes unbound, however, there is no controversy and 
the extrapolation errors are under control. 

We now consider the chiral extrapolation of the 3-body 
spectrum. This could be calculated using an EFT with 
explicit pions. Alternatively, the chiral extrapolation can 
be calculated using the contact EFT of Ref. |5j]. The 
inputs required are a s , at, and A* as functions of m^, 
which can be calculated using an EFT with pions. For 
the inverse scattering lengths l/o s and l/at, we take the 
central values of the error bands obtained from the chiral 
extrapolation in Ref. ^5|. The dependence of A* on m„ 
could be determined from the chiral extrapolation of the 
triton binding energy using an EFT with pions, but this 
has not yet been calculated. Like the inverse scattering 
lengths, A* should vary smoothly with m n . For small 
extrapolations of from its physical value, we can ap- 
proximate A* by a constant. We use the value A* = 189 
MeV for m,r = 138 MeV obtained by taking the triton 
binding energy as the 3-body input. In Fig. [2 we show 
the 3-body spectrum in the triton channel as a function of 
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TOtt. The crosses give the binding momenta k = (mB^) 1 / 2 
of the physical deuteron and triton, while the dashed 
lines give the thresholds for decay into a nucleon plus 
a deuteron (left curve) or a spin-singlet deuteron (right 
curve) in the large-a approximation. The circles indi- 
cate the triton ground state and excited state. In the 
region near m w w 175 MeV where the decay threshold 
comes closest to the 3-nucleon threshold K — 0, the ex- 
cited state of the triton appears in the spectrum. The 
existence of this very shallow excited state is a hint that 
the system is close to an infrared limit cycle. 

If, as in the case illustrated by Fig-El the critical values 
m Wt t and m Wia at which l/a t and l/a s go to zero satisfy 
TOtt^s < TOjr^t, then the deuteron and spin-singlet deuteron 
are both bound in the region to^s < m^ < to^. Since 
the decay threshold never extends all the way down to 
the 3-nucleon threshold k = 0, there cannot be an in- 
frared limit cycle. We argue that it is possible, however, 
to tune the system to the critical RG trajectory for an 
infrared limit cycle by adjusting the up and down quark 
masses m u and md- In the next-to-leading order chiral 
extrapolation of Ref. |l5j. only quark mass dependent 
operators proportional to m u + m ( i enter. The extrap- 
olation in m^ can be interpreted as an extrapolation in 
the sum m u + m c i, with m u — md held fixed. Chang- 
ing m u — md changes the degree of isospin-symmetry 
breaking. Since the spin-singlet and spin-triplet channels 
correspond to isospin-triplet and isospin-singlet, respec- 
tively, the scattering lengths a s and a t will respond dif- 
ferently to changes in m u — m<j. By tuning both m u and 
rrid, it should be possible to make l/at and l/d s vanish si- 
multaneously: m^i = m WlS . This point corresponds to a 
critical RG trajectory for an infrared limit cycle. At this 
critical point, the triton has infinitely many increasingly- 
shallow excitations with an accumulation point at the 
3-nucleon threshold. As the excitations get more and 
more shallow, the ratio of the binding energies of suc- 
cessive states rapidly approaches Aq « 515. Only for 
the deepest states like the triton will there be significant 
deviations. Now consider a renormalization group trans- 
formation that integrates out energies above some scale 
A 2 /m,N- As A is decreased, the deepest 3-nucleon bound 
states are removed from the spectrum, leaving only those 
for which the deviations from the asymptotic ratio Aq 
are negligible. Thus a limit cycle with a discrete-scaling- 
symmetry factor Ao is approached in the infrared limit 
A^O. 

The error bands in Ref. 0] do not exclude the critical 
values iTiir t and m^^ from satisfying m^j < "i^^, in 
which case the deuteron and spin-singlet deuteron are 
both unbound in the region ro^t < m^ < m^^. In this 
region, the decay threshold extends all the way down 
to the 3-nucleon threshold n = 0. At either of the two 
critical values and to^s, there is a 2-nucleon bound 
state at threshold. However there is no Efimov effect, 
because the 3-body problem is not sufficiently singular at 



short distances. The exact infrared limit cycle can only 
be obtained by using an additional tuning parameter such 
as m u — md to set m Vt t = m^.s- 

The infrared RG limit cycle of QCD has important 
implications for attempts to derive nuclear physics from 
lattice gauge theory. The computational effort for lat- 
tice simulations increases dramatically as the pion mass 
decreases and is prohibitive at the physical value. Lat- 
tice simulations are typically carried out at a value of ra, 
that is 2-3 times larger than the physical value, and then 
a chiral extrapolation is made to to m n = 138 MeV. Un- 
fortunately, this requires extrapolating past the region of 
m v where there is an RG limit cycle, which may intro- 
duce large extrapolation errors. The proximity of phys- 
ical QCD to the critical trajectory for the infrared limit 
cycle explains the success of Efimov's program 4] for de- 
scribing the 3-nucleon problem using zero-range forces 
or a contact EFT H . The apparent convergence of the 
effective-range corrections for momenta of order m, fl3| 
can be explained if the momentum expansion in the con- 
tact EFT is in powers of p/m*, where to* w 175 MeV is 
the critical value of the pion mass, which is significantly 
larger than the physical value. Our final conclusion is 
that the ultraviolet limit cycle in the contact EFT 
is not just an artifact of the model, but a hint of the 
existence of the infrared limit cycle of QCD. 
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